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1. INTRODUCTION 

Besides being a renewable and pollution free energy generation technology with no moving parts, 
PV modules can also be integrated into buildings as BIPV systems, adding aesthetic value [1]. When 
installed in an optimized way, BIPV systems can reduce heat transferred through the envelope and reduce 
cooling load components decreasing the CO2 emissions [2]. Apart from some facade installations, the 
rooftop segment represented more than 23 GWp of total installations in 2015, with projections of more than 
35 GWp to be installed by 2018 [3]. 

Since the BIPV offers the possibility to replace part of the traditional building material, with a 
possible price reduction in comparison to a classic rooftop installation [4], [5], the correct estimation of 
system level performances, system reliability and system availability is becoming more important and 
popular among installers, integrators, investors and owners; with this purpose several tools and models were 
developed [6-8]. The combination of different phenomena, such as the solar radiation available on site, the 
presence of dust, the shadowing or UV radiation over long outdoor exposure, affect in different ways the 
performance of BIPV systems and thus the related economic evaluations [9-11]. 

Many studies have been devoted to develop different non-linear electric models used to describe the 
characteristics of the PV modules and the effect on module performance of temperature, radiation intensity 
and other parameters and equipment/systems under non-standard conditions [12-20]. We offer a new method 
to model and analyze the performance of BIPV systems using Matlab/Simulink. 

In Section 2 of this article the mathematical description of the proposed model is presented to 
evaluate the performance of BIPV systems. Subsequently, Section 3 describes the 6 kW BIPV system 
installed. Section 4 presents the results obtained and the validation of the model with monitored data of the 
BIPV system. Finally, Section 5 presents the conclusions of the research carried out. 
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2. RESEARCH METHOD 
2.1. Photovoltaic Generator Model 
2.1.1. Solar Cell Model 
Figure 1 shows the model adopted for the solar cell [25-28]; which is widely used due to its 
simplicity and high degree of precision. 


Figure 1. Equivalent single-diode circuit of a solar cell 


Using circuit theory, the relation of the currents in the circuit can be obtained, which is given by: 
I = lpn — Ia — Ish (1) 


Where Iph represents the photo current generated which is directly proportional to the incident 
irradiance (G) and depends in less proportion on the operating temperature of the cell (T). This current can be 
expressed by the following equation: 


G 
Ton = Use + K; * (T = T,)] $ Gr (2) 


Where Isc is the short-circuit current under standard conditions (STC = 1000W/m”, 25°C, AM1.5) 
and Tr and Gr are the temperature and irradiance also in STC, Ki is the short-circuit temperature coefficient 
of the cell. 

The diode current Id represents the effect of the diffusion and recombination currents present at the 
PN junction of the cell. The most used mathematical model that allows to approximate the behavior of the 
diode is the Shockley equation that is given by: 


aVa 


la = Le * (elon) — 1) (3) 


In which q = 1.602x10~?° C is the charge of the electron, n is the diode ideal factor, dependent on 
the manufacturing process and usually adopts values between 1 (for germanium) and 2 (for silicon), 
k = 1.3806x107?3J/k is the Boltzmann constant, which relates absolute temperature and energy, Vd is the 
voltage at the terminals of the diode and Is is the saturation current of the cell, and is represented by the 
following equation: 


Ig = Ins * (2) el oe) re 


Where Epg = 1.12 eV is the bandwidth energy (in this case for silicon), and Irs is the inverse 
saturation current at standard temperature given by: 
Isc 


Ls = poza (5) 
e -1 


Ns*K*n>T. 


Where V,,. represents the cell’s open circuit voltage. 
Finally, Ish represents the leakage current through the junction and this is represented by the parallel 
resistance Rsh. 
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The leakage current is given by: 
ln (6) 
Where, 
Va =V +I* R; 0) 
Where V is the output voltage of the cell and Rs is the series resistance that represents the ohmic 


losses in the contacts due to the metal-semiconductor junction. Finally the following expression for Ish is 
obtained as: 


(8) 


Modifying expression (1) with Equations (3) and (8) gives the equation of the simple diode model 
for the solar cell to be used for the final model: 


V+Rs*I 
1 p=» (eF 1) — Pate (9) 


Rsh 


2.1.2. Photovoltaic Module Model 

The next step is to scale the model to an array because the power of a single cell is very small for 
general purposes. Generally the solar modules connect several cells in series (Ns) to produce more voltage 
and in parallel (Np) to produce more current. The following parameters were considered for scaling to a 
module: 


peni (10) 
Iscm = Ny * Isc (11) 
Va = Ne *V (12) 
Voom = Ns * Voc (13) 
Rsm = Rs = Ns (14) 
Rsnm = Rsn ba Ns (15) 


To represent a photovoltaic panel it is necessary to modify Equation (9) taking into account 
Equations (10) to (15) as follows: 


n*Vt*Ns 


N 
ae 


I = Ny * Inn — Ny * Ís * ( (16) 


Taking into account that for the model to develop Np = 1 the expression can be reduced as follows: 


RshM 


V+Rsm*I . 
I = lpn — lys (ele z 1) AKER (17) 


The Rs and Rsh values are not supplied by the manufacturer, but an approximate value of these 
values can be obtained by the following expressions: 


R = Voc—Vmpp (18) 


Smax Impp 
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-mee R, (19) 


Shmin ~~ Ise—Impp 

By means of iterative methods it is possible to obtain more approximate values of these resistances 
starting from the expressions (18) and (19). Some authors recommend using values as Rs less than 0.01Q and 
values of Rsh greater than 500 Q. 

The Newton-Raphson method is used because of the exponential nature of Equation (17), which 
must be solved, by using numerical methods such as the one mentioned above [29]. The output variables of 
the photovoltaic panel are produced with a vector from the open circuit voltage (Voc) and the power obtained 
(P). 

Once the equation of G vs Vmpp is obtained, it is implemented in the algorithm and the model is 
scaled to a complete generator taking into account the following relations: 


Ig = Npa * Iu (20) 
Vo = Nsg * Vu (21) 
Pg = Ig * Vo (22) 


Where Ic, Vc and Po represent the current, voltage and power of the solar generator respectively, 
Nyg represents the number of modules connected in parallel and N,grepresents the number of modules 
connected in series. 


2.2. Inverter Model 

The model proposed by Castañer [25] shown in Figure 2 is taken as reference. 

It is important to clarify that the inverters have incorporated the maximum power point (MPPT) 
tracker described above but the present model has been simplified since the generator model already delivers 
the values of Vmpp and Impp which are received through the resistors R1 and R2. The product of the voltage 
values in the resistors replicates the values of Pmpp given by the irradiance and the temperature in a certain 
period of time. In turn this product must be multiplied by the value of the efficiency of the inverter (nf) 
provided by the manufacturer. The controlled voltage source V1 must provide the sinusoidal amplitude 
values of the output current whereby the previously calculated power value must be divided by the voltage 
value provided by the electrical grid at Vpys due to being treated in power at DC as shown in the following 
relation: 


V; x]. anf 
Irus = Reppert (23) 
VredRMS 


Impp linv 


Ind Out1 


In2 


Figure 2. Equivalent circuit of the inverter model 


Since it is an efficient current, it is necessary to obtain the peak current values to obtain the 
amplitude values of the source V1 as shown in Equation (24): 
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Vmpp*ľImpp*nf* VZ _ I 
VredRMS p=p 


V = (24) 

The amplitude values obtained in V1 are used in the controlled current source linv by multiplying 
them by a sinusoidal signal (V2) which has the frequency and phase of the current provided by the electrical 
grid and a amplitude of 1V in order to obtain the corresponding waveform of the AC output current of the 
inverter. 


2.3. Electrical Grid Model 

The model used to simulate the power grid can be seen in Figure 3. 

Itis an RLC circuit whose main characteristic must be having a power factor greater than or equal to 
0.85 as established by the IEEE 929-2000 standard [30]. To meet such a requirement, an impedance value is 
first set to the appropriate angle set by Equation (25). 


linv Ired 


Figure 3. RLC circuit representing the local electrical grid 


0 = arccos(0.85) = 31.7° = 30° (25) 


Taking expression (26) into account, a value is established for R = 3Q and X =V3. The values of the 
inductive reactance (X,) and capacitive (Xç) are calculated considering the expressions (27) and (28) at a 
frequency of 60Hz. 


Tang = &*2 _ 8 (26) 
R 3 

X. = — 27 

Ces 2*0*f*C ( ) 

X, = 2nfL (28) 


An inductance value of L = 2H is established for the purpose of having a positive reactive power and 
the expression (29) is used to obtain the value of the capacitance (C) where w = 27fL [43]. 


C = —+—~ = 3.526 * 1076F (29) 


w2L—-wV3 


Once the development of all the components of the interconnected system has been completed, they 
are integrated into a single model (Figure 4), which will allow obtaining the behavioral results that represent 
the interaction of the different functional blocks. 
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Figure 4. Complete model of the photovoltaic system interconnected to the grid implemented in Simulink. 


3. SYSTEM DESCRIPTION 

The grid—connected BIPV system installed at the Universidad de Bogota Jorge Tadeo Lozano 
includes a PV array of 24 modules of poly crystalline silicon (Trina Solar TSM-PA05.08), each one of 250 
Wp and an inverter Sunny Boy 5000-US model of 5000 W. 

Taking into account that the DC input of the SB 5000-US inverter varies between 175 V and 480 V 
and the voltage at maximum power point (VMPP) of the module is 38V, the PV array was built 
interconnecting 2 branches in parallel of 12 modules in series each one. Under these conditions the nominal 
power of the PV array is 6000 Wp. 

The PV system is fully monitored to evaluate and analyze the energy produced and the power 
quality, since September 2015. The electrical variables are acquired every minute for every day and are 
stored in a spreadsheet file and they will be compared with the simulation results using Matlab. 


4. RESULTS AND ANALYSIS 
4.1. Photovoltaic Generator 

In order to evaluate the different characteristics of the module, the simulation is performed under 
different environmental conditions. Then, we maintain the operating temperature of the module (T) constant 
and we vary the irradiance (G) at different levels; the results obtained are shown in Figure 5(a) and 
Figure 5(b). 

As can be seen in Figure 5(a) and Figure 5(b), there is a strong relationship of the short-circuit 
current (Isc) with the increase of the incident irradiance (G). Open circuit voltage (Voc) is less sensitive. 


call 250; 

8l _200; 
<r =G = 
#6 Ta 
P j, 
ot a 

2: 


5 23 3 | 5 15 25 
i age (V) Voltage (V) 


(a) (b) 


Figure 5. (a) I-V curves, effect of irradiance at constant temperature of 25°C and (b) P-V curves, effect of 
irradiance at constant temperature of 25°C 
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For 800 W/m”, the current recorded was 7.2A with an open circuit voltage of 37.5V. For standard 
solar radiation of 1000 W/m’, the current obtained was 8.86A and an open circuit voltage of 37.8V. For the 
same level of solar radiation, the solar panel manufacturer reported a short-circuit current of 8.91A and open 
circuit voltage of 38.1V. This difference represents an error of 0.5% between the data generated by the model 
developed in Matlab™ and data reported by the manufacturer. 

Maximum power point in P-V curves: 

a. For a solar radiation level of 200 W/m? and 25°C of temperature, the open circuit voltage was 

29.5V and 48W of DC power. 

b. For a solar radiation level of 1000 W/m’ and 25°C of temperature, open circuit voltage was 30V 

and 248W of DC power. 

c. The percentage error between the maximum power point of the solar module reported by 

manufacturer and delivered by the model with 1000W/m’ and 25°C, was 0.8%. 

Figure 8(a) shows the comparison between the DC voltage of the photovoltaic generator measured 
by the monitoring system with the same voltage delivered by the model in Matlab. Figure 8(b) shows the 
comparison between the AC power of the photovoltaic generator measured by the monitoring system with 
the AC power delivered by the model in Matlab. 

The Matlab model has a voltage of 300V at 6:39 am; while the monitoring system registers a value 
of 355 V for the same hour. The maximum voltage reached by the solar panels is presented at 5:05 pm with a 
value of 410V. The correlation coefficient calculated between the DC voltage measured and that delivered by 
the model was 0.96. 


500) ___Measured DC Voltage —— Measured AC Power 
— Calculated DC Voltage 5000 Calculated AC Power 
400-4 
= | 
o 3004 
QD 
8 4 
2 200+ 
re) 
a 
100- 
0+ rere te ee eee ee Oo ttt ttt 
6 8" id 12 14 "ie ‘is 1 
Time (h) Time (h) 
(a) (b) 


Figure 8. (a) Monitored DC voltage and DC voltage of the model on a sunny day and (b) Monitored AC rms 
power and AC rms power of the model on a sunny day 


5. CONCLUSIONS 

In the present work the proposed model of a photovoltaic system interconnected to the electrical grid 
has been developed using the MATLAB/Simulink tool. The study has been carried out including the 
equivalent circuits of the fundamental components of its basic structure and analyzed under monitored 
irradiance conditions supplied by a measurement unit located in the facilities of an installed system. 

The complete model has been analyzed by means of a comparison of the nominal values of the 
modules and the inverter given by the manufacturers with respect to the data delivered by the model and a 
high degree of approximation was found, which suggests an adequate behavior compared to the standard 
conditions. 

The values of the model variables were compared with monitored measurements made at the 
Universidad de Bogota Jorge Tadeo Lozano with different irradiance characteristics and correlation 
coefficients with values higher than 0.92 were found, which confirms an adequate behavior of the model, 
very close to the monitored behavior of the pv system installed. 
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